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We obtain the entropy of flexible linear chains composed of M monomers placed on the square 
lattice using a transfer matrix approach. An excluded volume interaction is included by considering 
the chains to be self-and mutually avoiding, and a fraction p of the sites are occupied by monomers. 
We solve the problem exactly on stripes of increasing width m and then extrapolate our results 
to the two-dimensional limit m — » oo using finite-size scaling. The extrapolated results for several 
finite values of M and in the polymer limit M — + oo for the cases where all lattice sites are occupied 
(p = 1) and for the partially filled case p < 1 are compared with earlier results. These results are 
exact for dimers (M = 2) and full occupation (p — 1) and derived from series expansions, mean-field 
like approximations, and transfer matrix calculations for some other cases. For small values of M, 
as well as for the polymer limit M — > oo, rather precise estimates of the entropy are obtained. 



I. INTRODUCTION 

The term dimer was introduced in the thirties Q] as 
an abbreviation for diatomic molecules in a model for 
their adsorption on crystal surfaces. Later dimer models 
were applied in the study of other physical systems such 
as ferroelectrics, and much is known about their thermo- 
dynamic properties 2\. One relevant question in these 
models is the entropy associated to placing dimers on a 
regular lattice. For the particular case of full covering of 
the square lattice by dimers, this question was answered 
exactly some time ago, using a technique based on pfaf- 
fians [J. |j. y| However, even the generalization of this 
problem for the case on partial covering of the square lat- 
tice is still an open question today, no exact result being 
known. 

In this paper we address a generalization of the entropy 
of dimers problem, considering entropy related to cover- 
ing the square lattice with chains with M monomers each 
(we will call them M-mers) , as a function of the fraction p 
of sites of the lattice occupied by monomers. The chains 
will be considered flexible, so that there is no energy as- 
sociated to bending them. Since the only energy in the 
model is the infinite excluded volume interaction, which 
forbids the presence of more than one monomer on the 
same lattice site, the problem is athermal. It may be a 
simple model for the adsorption of monodisperse flexible 
chains on the surface of a crystal. Besides the exact so- 
lution of the problem for M = 2 and p = 1 mentioned 
above, other cases were already considered in the liter- 
ature. Rather precise transfer matrix calculations were 
performed in the polymer limit M — » 00 for hamiltonian 
walks (o = 1) @. There are also mean-field approxima- 
tions 0|, Bethe - and Husimi lattice calculations @, and 
series expansions in where q is the number of first 



neighbors of each site in the lattice Q , and in those calcu- 
lations approximate values for the entropy are obtained 
for both the full (p = 1) and partial (p < 1) coverage 
cases. 

In this paper we obtain estimates for the entropy of 
flexible M-mers placed on the square lattice, using trans- 
fer matrix techniques. This is done calculating numeri- 
cally exact values for the entropy of the problem on strips 
with finite widths m and periodic boundary conditions 
and then using finite-size scaling to extrapolate the re- 
sults to the two-dimensional limit m — > 00. We separate 
the problem in cases where M is finite or infinite (poly- 
mer limit). Also, the case of full coverage (p = 1) may be 
treated separately from the general case. In the general 
case, it is convenient to address the problem in an ensem- 
ble which is grand-canonical with respect to the number 
of monomers placed on the lattice, whereas for full cover- 
age it is easier to perform a microcanonical calculation. 

The expressions we used to calculate the entropy are 
shown in section^ The model is discussed in more detail 
and the transfer matrices are described in section IIIII 
Our results for the entropies may be found in section IWl 
as well as the extrapolation procedure and their results. 
Section presents final discussions and conclusions. 



II. DETERMINATION OF THE ENTROPY 

For the case of full coverage, it is convenient to obtain 
the entropy directly from Boltzmann's expression 
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where fl is the number of ways to fill the lattice with 
N sites completely with M-mers. In the polymer limit 
M — > 00, we consider a single hamiltonian walk, that is, 
a self-avoiding walk (SAW) which visits all sites of the 
lattice. 

In the general case where a fraction p of lattice sites 
is occupied by monomers, we define the grand-canonical 
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partition function 

E(z) = ^2^ M r{M,N,p), 



(2) 



where z is the activity of a monomer and T(M, N, p) is the 
number of ways to place p chains with M monomers each 
on the lattice with N sites. For the polymer limit, again 
a single SAW is placed on the lattice and the partition 
function is defined as 



S(*) = 5>T(n,i\0, 



(3) 



where T(n, N) is the number of ways to place a SAW 
with n monomers on the TV-site lattice. The density of 
monomers may now be written as 



(4) 



where the thermodynamic potential per lattice site is de- 
fined as 



(z) = lim — lnE(z). 

N— >oo 1\ 
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In the thermodynamic limit, a Legendre transformation 
allows us to rewrite the potential as 



4>{z) ~ maxjplnz + s(p)}, 
p 

and thus the entropy may be written as 

s(p) = - [ In z(p')dp', 
Jo 

with s(0) = 0. 
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2. The pairs of bonds which are connected to each 
other through a path lying entirely below the refer- 
ence line (see figure QJ. These pairs may be also be 
specified by a m-component vector |i>), associating 
a different positive integer to each pair of connected 
bonds and to the ones which are not connected 
to any other. This connectivity information pre- 
vents us from closing a ring at any level, since this 
configuration is not allowed in the model. 
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Figure 1: Example of a state for m = 4 vertical bonds. The 
reference line indicates the set of vertical bonds whose config- 
uration is described. 



As an example of these definitions, we consider the case 
of pentamers (M — 5) placed on a strip of width m = 4. 
Among the possible configurations of a set of vertical 
bonds, the one depicted in figure is described by the 
vectors \v) = (0,1,1,0) and \p) = (1,2,2,0). Elements 
of the line associated to this state of the transfer matrix 
T are obtained considering the possible continuations of 
the state {\v ), \p)} one step upwards, as shown in two 
examples in figure El The resulting state in figure EJa 
is described by the vectors \v) = (0,0,0,0) and \p) = 
(2,0,4,0), while the final state in figureEJb corresponds 
to the vectors \v) = (0, 0, 0, 0) and \p) = (2, 0, 0, 0). Each 



III. DEFINITION OF THE TRANSFER 
MATRIX 

We proceed defining a strip of width m on the square 
lattice in the (x, y) plane, so that 1 < x < m and 
— oo < y < oo, with periodic boundary conditions in both 
directions. A transfer matrix may be built for this prob- 
lem, inspired on the prescription due to Derrida [lfj] for 
infinite chains in strips. We thus consider the operation 
of including an additional step to the strip in the positive 
y direction, adding m new sites of the lattice. To properly 
take into account the statistical weight of the new step, 
we may define the state of the m vertical bonds of the 
lattice which are incident to the new sites by specifying: 

1. The number of monomers already present in the 
chain which passes through the vertical bond (it is 
equal to if no chain is present). These numbers 
may be put into a vector \p), with m components. 
It is necessary to keep track of this information so 
that we know when to end each chain. 




(a) (b) 

Figure 2: Possible continuations (reference line RL') following 
the configuration depicted in U (reference line RL). 

monomer placed on a site between the two sets of ver- 
tical bonds contributes with an activity z to the parti- 
tion function, so that the element of the transfer matrix 
which corresponds to the first configuration is equal to 
z 3 , while the second configuration is associated to an ele- 
ment equal to z A in the transfer matrix. Only the second 
configuration contributes in the case of full occupancy. 

For the polymer case M — > oo, a single chain passes 
through the whole strip, so that it is enough to describe 
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the connectivity at a particular set of m vertical bonds 
by indicating the bond which is connected to the initial 
monomer of the chain (in y — » — oo) and the pairs of 
bonds connected to each other, exactly as was done in 
the original work of Derrida [1(T| . Thus, a single vector 
\v) is enough to describe the state in this limit. 

Once the transfer matrix T is obtained, the entropy of 
the model on the strip in the thermodynamic limit is re- 
lated to the largest eigenvalue of the matrix. For the case 
of full occupancy (p = 1), the number of configurations 
is given by 



n = Tr(T' 1 ), 



(8) 



where N = ml is the number of sites and the elements of 
the matrix T' are defined by the limit 
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The entropy is then related to the largest eigenvalue A 
of this matrix, so that 
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For the general case, where a fraction p of lattice sites 
are occupied by monomers, the grand-canonical partition 
function is related to the transfer matrix through 



S(z)=Tr(T l ), 
and thus the density p(z) will be 

p(z) = lim ^-^-lnE[z) = 

n-kx> N dz 
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where A is the largest eigenvalue of the transfer matrix 
T. This relation may be inverted to obtain the entropy 
as a function of the density using equation 



IV. NUMERICAL RESULTS 

The size of the transfer matrix increases very fast with 
both the molecular weight M and the width m of the 
strip grow. This sets an upper limit to the widths we 
were able to consider for each chain of a given molecular 
weight. Figure El shows this effect, and in the inset one 
may appreciate that the growth of the transfer matrix is 
roughly exponential. 

Furthermore, as was already observed in similar cal- 
culations for polymers the values of the entropy for 
each class of M-mers are split into subsets with differ- 
ent finite-size scaling behavior in each subset according 
to the width of the strips, so that extrapolations must 
be done within each subset. These splitting seems to be 
related to frustration effects in the limit of the fully oc- 
cupied lattice and the subsets are indicated in table Q] 
where all the widths we considered are given. 
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Figure 3: Dimension of the transfer matrix as a function of 
the molecular weight (m = 2) and width of the strip (M = 2) , 
respectively. 



The data of each subset were extrapolated to the two- 
dimensional limit to — > oo using the Shanks transfor- 
mation 0|, since we expect finite-size corrections to be 
exponential. Since at least three values for the entropy 
are needed in a subset to obtain an estimate for the two- 
dimensional entropy and its confidence interval, not all 
subsets may be extrapolated, so that, for example, no es- 
timate could be found for hexamers (M = 6), where we 
calculated entropies for widths up to m — 7. The final 
estimate was chosen to be the highest possible extrap- 
olant, and the error associated to it was obtained from 
the previous generation of extrapolants, through 



e = lim 
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The extrapolated values of the entropies for p = 1 and 
their uncertainties are displayed in table[nj together with 
values obtained with other techniques and best values 
found in the literature. 

Our results may be compared with other values in the 
literature. One may notice that the mean-field estimates 
are systematically smaller than the values obtained here, 
but no such general trend is apparent for the Bethe- and 
Husimi lattice results. The estimate for dimers agrees 
with the exact value obtained by Fisher, Temperley and 
Kasteleyn 0,0,0], and the entropy for hamiltonian walks 
(M — > oo) is consistent with both transfer matrix calcu- 
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Table I: Entropies calculated for each M-mer, divided in sub- 
sets with the same finite-size scaling behavior 

Molecular weight Entropies of each subset 



M = 2 



{si} = {s 2 , S 4 , S 6 , S14} 
{s' 2 } = {S3,S5,S7, ...,S13} 



M = 3 



{si} = {S3, S6, S9, S12} 
{ s 2} = { s 2, 54, S5, S7, S8, SlQ, Sn} 



M = 4 



{si} = {s 4 ,s s } 

{s 2 } = {S2, S6, S10} 
{ s 3 } = {S3,S 5 ,S 7 ,S 9 } 



M = 5 



{si} = {«} 

{s 2 } = {S2, S3, S4, S6, S7, Ss} 



M = 6 


{si} = {s 6 } 

W2} = {S3} 
{S3} = {S 2 ,S 4 } 
{s 4 } = {S 5 ,S 7 } 




M = 7 








{si} = {s 2 , S3, S 4 , S 5 


s 6 } 


M = 8 


{si} = {s 4 } 

{S' 2 } = {S2} 

{s 3 } = {s 3 ,s 5 } 




M = 9 


{Si} = {s 3 } 
{S 2 } = {S2,S 4 } 




M — > 00 








{si} = {S 2 ,S 4 ,S 6 , 


Sl 2 } 




{s 2 } = {s 3 ,s 5 ,s 7 , 


S13} 



lations |6| and the result of series expansions up to third 
order in q~ x 0], which is Soo w 0.38629. 

Another relevant question is the value of the molecular 
weight which maximizes the entropy at full occupancy. 
Mean-field and Bethe lattice results show maximum en- 
tropy at M = 8 for a lattice with coordination number 
9 = 4 JE 111; while series up to second order in g -1 on the 
square lattice result in a maximum entropy at M = 7. 
Our results suggest that this maximum actually occurs 
at M = 4 on the square lattice, if we suppose that that 
only one maximum exists in the curve s(p = 1) x M and 
also disregard the value obtained for M — 7 due to the 
large uncertainty associated to it. 

For partial occupancy of the lattice, the results are 
similar to the ones shown in figure for the entropy of 
dimers as a function of the fraction of occupied lattice 
sites p. For all the cases we considered the entropy dis- 
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Figure 4: Entropy at full occupancy of the lattice as a function 
of the molecular weight M 



plays a single maximum. The density at which this maxi- 
mum occurs increases with M, getting closer to p « 0.79, 
the value found in the polymer limit. The densities and 
maximum entropies are listed in table Hill and it may be 
noticed that the largest value for the maximum entropy 
occurs for tetramers, as was also found for p = 1. 
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Figure 5: Entropy for dimers as a function of the density. In 
the inset, results of values obtained in this work are compared 
with results from series expansion Q and Bethe lattice [j| 
calculations. 



In the polymer limit M — > 00, the model of a poly- 
mer placed on a strip in the grand-canonical ensemble 
displays a first order phase transition at a critical activ- 
ity z c , with the coexistence of a non-polymerized phase 
(p = 0) and a polymerized phase (p = p c > 0) [Tj|. As 
the width m of the strip is increased, the discontinuity in 
the density at the transition becomes smaller and in the 
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Table II: Entropy of M-mers on the square lattice for full coverage (p = 1) obtained through different techniques 



M 


MF 


Bethe 


Series 


Transfer Matrix 


Best value 


2 


0.19315 


0.26162 


0.26867 


0.29120± 0.00071 


0.29156 


3 


0.39268 


0.42284 


0.41699 


0.41201± 0.00002 




4 


0.46301 


0.48166 


0.48889 


0.51486± 0.0045 




5 


0.49229 


0.50669 


0.51008 


0.49917± 0.00091 




7 


0.51008 


0.52217 


0.52170 


0.54770± 0.15301 




M — > oo 


0.3863 


0.4055 


0.4090 


0.3870± 0.0009 


0.3866 



Table III: Maximum values of the entropy as a function of the 
density. 

M density of maximum entropy maximum entropy 



2 
3 
4 
5 
7 

M — > oo 



0.64 
0.71 
0.76 
0.76 
0.78 
0.79 



0.66 
0.70 
0.74 
0.73 
0.72 
0.56 



two-dimensional limit to — > oo a continuous transition 
is found at z c w 0.3790522 The entropy of a poly- 
mer on a strip of finite width is therefore not defined for 
p < p c , as may be seen in figure® where the extrapolated 
value of the entropy as a function of the density for poly- 
mers is depicted. As larger widths are considered, the 
step in the entropy decreases and in the two-dimensional 
limit we have s(p = 0) =0. 



V. CONCLUSION 

In this paper we estimate the entropy of chains with M 
monomers each placed on the square lattice as a function 
of the fraction p of lattice sites occupied by monomers. 
The estimates were obtained by extrapolating numeri- 
cally exact values for the entropy on strips of finite widths 
to, calculated using a transfer matrix approach, to the 
two-dimensional limit to — > oo. 

A particular case where results may easily be obtained 
analytically is the one dimensional problem (m = 1). In 
this case, the dimension of the transfer matrix T is equal 
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Figure 6: Entropy of polymers as a function of the density. 
The squares correspond to extrapolations of the values on 
strips of finite widths and the full line is the result for a strip 
of width m = 7. 



to M, and we have 

%,j = Si,i6j,i + z(8i+\j + 5i,MSj,i), (14) 

where 1 < i,j < M. One may easily obtain the secular 
equation of this matrix, which is 



A M _ A M-1 _ Z M = _ 



(15) 



The density may then be found as a function of the 
largest eigenvalue A through using the secular equation 
above and also expression El One obtains 



A 



1 — ap 



(16) 
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where a = (M — 1)/M. Then the entropy may be found 
by performing the integration in equation changing the 
integration variable from p to A. The result is 

s = (1 - ap) ln(l - ap) - (1 - p) ln(l - p) + 

- p(l-a)Hp(l-a)]. (17) 

This result is equal to the expression which is obtained 
if the coordination number of the Bethe lattice result 
(expression 22 in @) is taken equal to 2. The entropy in 
the one-dimensional case vanishes, as expected, for p = 1, 
and the maximum is located at a value of the density 
which is equal to 1/2 for monomers (a = 0) increasing 
monotonically with a and approaching 1 in the polymer 
limit a — ► 1 (figure EJ| , where the entropy vanishes for 
all values of p. The value of the maximum entropy is a 
decreasing function of a. 
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Figure 7: Density of maximum entropy as a function of 1/M. 
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For a given value of M, the density which maximizes 
the entropy is obtained through he equation 



(l-ap)- a (l-pMl-a)}-^ = l 



(18) 



On the square lattice, our calculations indicate that 
the absolute maximum of the entropy sm(p) occurs for 
tetramers (M = 4). In the one-dimensional case, the 
maximum entropy of sm{p) is a monotonically decreas- 
ing function of M, the absolute maximum si(l/2) = 
ln(2) being obtained for monomers and approaching 
as M — > oo. 

Finally, the problem of M-mers confined inside strips 
of finite width with closed boundary conditions, is an 
interesting extension of earlier work done in the polymer 
limit 0. 
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